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Abstract 

Bobecka and Wesolowski (2002) have shown that, in the Olkin and Rubin characterization of 
the Wishart distribution (See Casalis and Letac (1996)), when we use the division algorithm defined 
by the quadratic representation and replace the property of invariance by the existence of twice 
differentiable densities, we still have a characterization of the Wishart distribution. In the present 
work, we show that, when we use the division algorithm defined by the Cholesky decomposition, 
we get a characterization of the Riesz distribution. 

Keywords: Symmetric cone, division algorithm, Wishart distribution, Riesz distribution, Beta- 
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1 Introduction 

A remarkable characterization of the gamma distribution, due to Luckacs (1955), says that if U and 
V are two independent non Dirac and non negative random variables such that U+V is a.s. positive, 
then jj^y and U + V are independent if and only if U and V have the gamma distribution with 
the same scale parameter. The classical multivariate version of this characterization concerns the 
Wishart distribution on the cone of symmetric positive matrices. In this case, there is not a single 
way to define the quotient of two matrices. For instance if Y is a positive definite matrix, one 
can for example, use the quadratic representation, that is write Y = Y2Y2 and define the ratio 
X by Y as Y~?XY~2, or use the Cholesky decomposition Y = TT* , where T is a lower triangular 
matrix and define the ratio as (T~ 1 )X(T~ 1 )* . A general definition of a division algorithm will 
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be given in Section 2, however these two examples are the most usual and most important. In 
1962, Olkin and Rubin have shown that, if U and V are two independent random variables valued 
in the cone of symmetric non negative matrices such that U + V is a.s. positive definite, then, 
independently of the choice of the division algorithm, the quotient of U by U + V is independent 
of U + V and its distribution is invariant by the orthogonal group if and only if U and V have the 
Wishart distribution. This result has been extended to the Wishart distribution on any symmetric 
cone by Casalis and Letac (1996). Recently, Bobecka and Wesolowski (2002) have given another 
characterization of the Wishart distribution without any invariance assumption for the quotient. 
More precisely, they have shown that if we use the division algorithm defined by the quadratic 
representation and we replace in the Olkin and Rubin theorem the condition of invariance for 
the distribution of the quotient by the existence of twice differentiable densities, then we still 
have a characterization of the Wishart distribution. The present paper gives a parallel result 
which starts from the observation that, when the condition of invariance of the distribution of 
the ratio by the orthogonal group is dropped, the characterization in the Bobecka and Wesolowski 
way is not independent of the choice of the division algorithm. We show that, when we use the 
division algorithm defined by the Cholesky decomposition, we get a characterization of the Riesz 
distribution introduced by Hassairi and Lajmi (2001). Our method of proof is based on some 
functional equations depending on the triangular group. These equations are more involved then 
the ones used in the characterization of the Wishart distribution, their solutions are expressed 
in terms of the generalized power. Our results will be presented in the framework of the Riesz 
distribution on the symmetric cone of a simple Euclidean Jordan algebra. This will enable us to 
use some technical results established in the book of Faraut-Koranyi (1994) and in Hassairi et al. 
(2001, 2005). However, to make the paper accessible to a reader who is not familiar with the theory 
of Jordan algebras, we will give a particular emphasis to the cone of positive definite symmetric 
matrices. 

2 Riesz distributions 

We first review some facts concerning Jordan algebras and their symmetric cones. Our notations 
are the ones used in the book of Faraut-Koranyi (1994). Let us recall that a Euclidean Jordan 
algebra is a Euclidean space E with scalar product < x, y > and a bilinear map 

E x E — > E, (x,y)\ — ► xy 

called Jordan product such that, for all x, y, z in E, 

i) xy = yx, 

ii) < x,yz > = < xy, z >, 

iii) there exists e in E such that ex = x, 

iv) x(x 2 y) = x 2 (xy), where we used the abbreviation x 2 = xx. 

An Euclidean Jordan algebra is said to be simple if it does not contain a nontrivial ideal. 
Actually to each Euclidean simple Jordan algebra, one attaches the set of Jordan squares 

H= {x 2 ; x e E}. 
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Its interior f2 is a symmetric cone, i.e., a cone which is 

i) self dual, i.e., = {x £ E; < x,y > > Vy eQ\ {0}}. 

ii) homogeneous, i.e., the subgroup G(Q) of the linear group GL(E) of linear automorphisms which 
preserves £1 acts transitively on Q. 

iii) salient, i.e., f2 does not contain a line. Furthermore, it is irreducible in the sense that it is not 
the product of two cones. 

Let now x be in E. If L(x) is the endomorphism of E; y i — ► xy and P(x) = 2L(x) 2 — L(x 2 ), 
then L(x) and P (x) are symmetric for the Euclidean structure of E and the map x i — > P(x) is 
called the quadratic representation of E. 

An element c of E is said to be idempotent if c 2 = c, it is a primitive idempotent if furthermore 
and is not the sum i + u of two non null idempotents t and u such that t.u = 0. 

r 

A Jordan frame is a set {c\ , c r } of primitive idempotents such that ^ a = e and CjCj = <%q, 

i=l 

for 1 < i,j < r. It is an important result that the size r of such a frame is a constant called the 
rank of E. 

If c is a primitive idempotent of E, the only possible eigenvalues of L(c) are 0, \ and 1. The 
corresponding eigenspaces are respectively denoted by E(c,0), E(c, ^) and E(c, 1) and the decom- 
position 

E = E(c,0) ® E(c, -) ® E(c,l) 

is called the Peirce decomposition of E with respect to c. 

Suppose now that {ci)i<i<r is a Jordan frame in E and let, for 1 < i,j < r, 

f E(ci,l) = m* if i = j 

^ \ £(Q,±)n£(c„±) if^j. 

Then (See Faraut-Koranyi (1994), Theorem IV.2.1) we have E = ®i<jEij and the dimension of Eij 
is, for i ^ j, a constant d called the Jordan constant. It is related to the dimension n and the rank 
r of E by the relation n = r + r(r — 1)|. 

For 1 < k < r, let denote the orthogonal projection on the Jordan subalgebra 

EW = E( Cl + ... + c k ,l), 

det^-* the determinant in the subalgebra E^ and, for x in E, A fc (x) = det^(P fe (x)). Then 
Afc is called the principal minor of order with respect to the Jordan frame (cj)i<j< r . For s = 
(si, s r ) £ IR r , and x in 17, we write 

A s (x) = A!(x) Sl - S2 A 2 (x) S2 - S3 A r (x) a -. 
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This is the generalized power function. Note that, if x = AjCj, then A s (x) = A* 1 AJ 2 ■ ■■Xp' and that 

i=i 

A s (x) = (detx) p if s = (p, ...,p) with p£E It is also easy to see that A s+s /(x) = A s (x)A s r(x). 
In particular, if m £ 1R and s + m = (si + m, s r + m), we have A s+m (x) = A s (x)(det x) m . 

As we have mentioned above, one may suppose that E is the algebra of real symmetric matrices 
with rank r. In this case, the Jordan product xy of two symmetric matrices x and y is defined 
by \{x.y + y.x) where x.y is the ordinary product of the matrices x and y, the cone 17 is the 
cone of positive definite matrices, 17 is the cone of symmetric non negative matrices and d = 1. If 
x = {xij)i<i,j< r is an (r, r)— symmetric positive definite matrix, and if, for 1 < k < r, we denote 
Pk{x) = (xij)i<i j<k and Afc(x) = det(xjj)i<j j<&, the generalized power is the function on f2 defined 
by 

A s (x) = A 1 (x) s i- s *A 2 (x) s *- S3 A r (x) s r. 



The definition of the Riesz distribution on a symmetric cone 0, relies on the notion of generalized 
power. In fact, for a is in f2 and s such that, for all i, Si > (i — 1)|, the measure on 

R(s,a)(dx) = ——^-—e- < ^ > A s ^(x)l Q (x)dx 

r ^ 

where Tq(s) = (2ir)^~ T(sj — (j — 1)-), is a probability distribution. It is called the Riesz 

3=1 2 
distribution with parameters s and a. 

We come now to the general definition of a division algorithm in a symmetric cone fL Let G be 

the connected component of the identity in G(Q). A division algorithm is defined as a measurable 

map g from 17 into G such that, for all y in 17, g(y)(y) = e. 



As in the case of symmetric matrices, we will introduce two important division algorithms, the 
first is based on the quadratic representation x 1— > P{x~z ) and the second algorithm takes its values 
in the triangular group T. For the definition of T, we need to introduce some other facts concerning 
a Jordan algebra. For x and y in E, let xOy denote the endomorphism of E defined by 

xUy = L(xy) + [L(x),L(y)} = L(xy) + L(x)L(y) - L(y)L(x). (1) 

If c is an idempotent and if z is an element of E(c, |), 

t c (z) = exp(2zDc) 

is called a Frobenius transformation, it is an clement of the group G. 

Given a Jordan frame (cj)i<j< r , the subgroup of G 

T=lr Cl (z^) r Cr _ 1 (^- 1 ))P(^a iQ ), a t >0, e ® E jk 

I i=l fc=j+l 
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is called the triangular group corresponding to the Jordan frame (cj)i<j< r . It is an important result 
(Faraut-Koranyi, p. 113, Prop VI. 3. 8) that the symmetric cone £1 of the algebra E is parameterized 
by the set 



E+ = {u = Y^ UiCi + u ij, u i > °} ( 2 ) 

i=l i<j 



More precisely, if 



t u = r ci (z {1) ) r Cr _ 1 (^- 1 ))P(^n,c i ) (3) 

i=i 

r 

where Zij = ^f, i < j and = ^ Zjk, then the map u i — ► i«(e) is a bijection from into 

k=j+l 

r 

with a Jacobian equal to 2 r TT u^^ 1 ^ . Also, for all x in £7, we have 

i=\ 

A fc (t„(s)) = u?...«gA fc (x) = A fc (* u (e))A fc (x). 

It is shown that, for each 6 in Q, there exists a unique t in the triangular group T such that 6 = t(e). 
Hence the map 

g : n — ► T; 6 i — ► i _1 (4) 

realizes a division algorithm. 

3 Characterization of the Riesz distribution 

In this section we state and prove our main result which may be seen as an extension of the result 
by Bobecka and Wesolowski (2002) concerning the ordinary Wishart distribution on symmetric 
matrices. More precisely these authors use the division algorithm defined by the quadratic rep- 
resentation to characterize the Wishart distribution, we use the division algorithm defined by the 
triangular group to characterize the Riesz distribution. 

Theorem 3.1. Let X and Y be two independent Riesz random variables X ~ R(s,a) and 
Y ~ R(s', a). If we setV = X + Y and U = g(X + Y)(X), then 

i) V is a Riesz random variable V ~ R(s + s', a) and is independent of U 

ii) The density of U with respect to the Lebesgue measure is 

where Bq(s,s') is the beta function defined on the symmetric cone Q (See Faraut-Koranyi, 1994, 
p. 130) by 

r n ( s )r n (s') 
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Proof. Consider the transformation: S7 x $7 — ► (17 n (e — J7)) x $7; (x, y) i — ► (u,v), where 
v = x + y = te, t G T and u = Its Jacobian is deti -1 = (dett(e))~~ = (detv)~~ . 



The density of probability of (U, V) with respect to the Lebesgue measure is then given by 

fgi^P) A.-;W»))A,-- ; (t(e-»)) (det^ «-<^U(«,,) 
where K is defined by 

K = {(u, v)/u £ Qn(e- Q) and v G fi}. 
Using equality (3.5) in Hassairi and Lajmi (2001), this density may be written as 

Jrc/fJ" ^ A s _ 2 (u) A s ,_n(e-u)— _i — _ e ~ <a ' v> A s+s ,_n(v)l K (u,v). 

From this we deduce that U and V are independent and that V ~ i?(s + s',a). Furthermore the 
distribution of U is concentrated on Q n (e — $7) with a density equal to 

Tq(s + s') 

□ 

Note that the distribution of the random variable U = g(X + Y)(X) is called beta-Riesz 
distribution with parameters s and s' (See Hassairi et al (2005)). In the case of symmetric matrices, 
the random variable U is nothing but 

U = (T-^XiT- 1 )*, 
where T is a lower triangular matrix with positive diagonal such that 

X + Y = TT*. 

Theorem 3.2. Let b i — ► gib) be the division algorithm defined by (4). Let X and Y be independent 
random variables valued in S7 with strictly positive twice differentiable densities. Set V = X + Y and 
U = g(V)(X). IfU and V are independent then there exist s, s' G IR r ; Sj > (i— s[ > (i — 1)| 
for all i, and a G S7 such that X ~ R(s, a) and Y ~ R(s', a). 

The proof of this theorem relies on the resolution of two functional equations given in the fol- 
lowing theorems which are interesting in their own rights. The proofs of these theorems are given 
in Section 4. 

Theorem 3.3. Let a : £1 n (e — J7) — ► IR and g : Q — ► IR be functions such that, for any 
x G Q n (e - £7) and t G T, 

a{x) = g(tx) - g(t(e- x)). (5) 
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Assume that g is differentiable, then there exist p £ IR r and c £ IR such that, for any x G Hfl(e- fi) 
and y £ 

a(x) = log A p (x) - log A p (e - x), g(y) = log A p (y) + c. 



Theorem 3.4. Zei a\ : £1 l~l (e — 17) — ► IR and 02, 3 : fi — ► IR be functions satisfying 

ai(x) + a 2 (te) = g(tx) + g(t(e - x)), (6) 

/or any x £ Qf] (e — Q) and t £ T . Assume that g is twice differentiable then there exist p' £ IR r , 
5 £ E and c±, c 2 , C3 £ IR such that for any x £ Q, n (e — 0) and y G fl, 

3(y) = logA p /(n)+ < (5,2/ > +ci 

ai (x) = log A p / (x) + log A p i (e — x) + c 2 , 
a 2 (y) = 21ogA p /(n)+ < <5,y > +c 3 , 

where 2c\ = 02 + 03. 

Proof of Theorem 3.2. We again use the transformation: x Q — ► (fin (e — Q)) x fi; (x, y) 1 — ► 
(u, t>), where v = x + y = te, t £ T and n = t~ 1 (x). Let /x, /y, fu an d /y be the densities of 
X, Y, U and V, respectively. Then, since (X, Y) and (U, V) have independent components, we 
have that, for all u £ £1 n (e — f2) and 

fu(u)fv(v) = (dett;)? fx(tu)f Y (t(e - «)) (7) 
Taking logarithms in (7) we get 

yiO) +52(te) = 33^) +S4(*(e - it)), (8) 

where 



31 (u) = log/c/(w), 



n 



32(f) = log/y(u) logdetu, 



33 (x) = logfx(x), 



(9) 
(10) 

(11) 



and 



34(3) = log My). (12) 

Inserting e — n for n in (8) gives 

31 (e -u) + g 2 {te) = g^{t(e - u)) + 34 (iu). (13) 
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Subtracting (13) from (8), we obtain 

9i H -gi(e-u) = [gs(tu) - g A {tu)) - \g 3 (t(e - u)) - g A {t{e 

Define 

a ( u ) = gi( u ) - gi( e - u ), 9 = 53 - 54, 

then, 

a(u) =g(t(u)) -g(t(e-u)). 
Now according to Theorem 3.3, we obtain 

a(u) = log A p (u) - log A p (e - u), g(v) = log A p (v) + c 

for p = (pi, ...,p r ) G lR r and c € IR. 
Hence 

53 (w) = 54 (v) +g(v) = g±(y) + logA p (v) + c. 
Inserting (15) back into (8) gives 

gi(u) + g 2 (te) = logAp(u) + logA p (ie) + c + g A {t{u)) + g A {t{ 

which can be rewritten in the form 

ai(u) + a 2 (te) = g A (t(u)) + g A (t(e - u)), 

where 

ai(u) = gi{u) - log A p (u), 

a 2 (te) = g 2 (te) - log A p (te) - c. 
Hence, by Theorem 3.4, it follows that 

g A ( V ) = log A p /(t;)+ <S,V>+Ci, 

ai(u) = log Apt (u) + log Apt (e — u) + c 2 , 
a 2 (-u) = 21ogA p /(i;)+ < S,v > +c 3 , 
for some 5 e E, p' <E IR r and c\, c 2 , C3 G IR such that c 3 = 2c\ — c 2 . 

This with (12) imply 

log /y(y) = logA p /(y)+ <5,y>+c 1 , 

that is 

/y(5) = e Cl e<^>A p ,(5). 



Since fy is a probability density, it follows that a = —5 £ £1 and s' = p' + j is such that 
s'i > (i - 1)| , VI < i < r so that F ~ R(s', a). 



Now by (15), we get 



g 3 (v) = logA p (f) +c + logA p /(w)+ < > +c x 
= log A p+J) /('y)+ < <5, v > +c + c\. 

From (11) it follows that 

f x (x) = A p+P ,(x)e«^> e c+Cl 
which implies that pi + s'i > (i — 1)|, VI < i < r and consequently X ~ i?(s, a) where s = p + s'. 

4 Proofs of Theorems 3.3 and 3.4 
4.1 Proof of Theorem 3.3 

For the proof of Theorem 3.3, we need to establish some preliminary results. 
Proposition 4.1. i) The map (p : £1 — > H; x i — ► log A^(x) is differentiable on and 

VlogA k (x) = (P k (x))- 1 . 
ii) The differential of the map x i — ► (Pk{x))~ l is —P((Pk{x))~ l ). 
Proof. We first observe that if c is an idempotent of E and x is in E(c, 1), then 

P{x)P{c) = P{x) 

In fact, let h = h\ + h\i + ho be the Peirce decomposition of an element h in E with respect to c. 
As x e E(c, 1), we have that 

P(x)(h 12 ) = P(x)(h ) = 0, 

and it follows that 

P(x)(h) = P0r)(M = P(x)P(c)(h) 

i) Let Qk = -Pfc(^) and consider the map ip : VL^ — ► IR; x \ — ► logdet^(x). Then tp(x) = tpoP k (x). 
Since Vlogdetx = we have 
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(logA k (x)Y(h) = ^(P k (x))(P k (x))\h),VheE 
= <(P k (x)Y\P k (h)> 
= <(P k (x))-\h> 

ii) We have that the differential in x of the map (3 : x i — ► x _1 is — P(x _1 ). As (P k (x))~ 1 = (3 o 
Pfe(x), Vx G 0, then the differential in x G f2 of the map a; i— > (i-^x)) -1 is equal to — P((Pfc(x)) _1 )o 
Pfc which is also equal to — P((P k (x))~ 1 ) 

Proposition 4.2. Lei u be in E + and let x = t u (e). Then, for 1 < k < r, 

k 

{p k (x)y l = ti\Y.^). 

i=l 

Proof. We know, from Faraut-Koranyi (Proposition VI. 3. 10), that for x in E and t u in T, 

P k {tu{x)) = t Pk{u) {P k {x)). 

Using (3) we can write 

P k (x) = T(z^)...r(z^)P(j2u l c t )(c l +... + c k ) 

i=i 

i=i 

k 

where = z jt and z jt = ^ L , j < I- 
l=j+i 

Hence 



(Pkix))- 1 = r(-z^y...T{-z^T (£ \«) (19) 

i=i u i 

On the other hand 

t* u l ctc l ) = T{-z^y...T{- z ^)y (£±1*). 

i=l i=l U i 

Let us recall that if c is an idempotent, z is in E(c, ^) and x\, xu, %o are the Peirce components 
of x with respect to c. Then the Peirce components of y = t c (z)*(x) are 

yi = 2c[z(zx ) + ZX12] + xi, 

yi2 = 2zx + X12, (20) 
2/0 = zo- 
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This, after some elementary calculations, implies that, for k + 1 < j < r — 1, E E(cj, |) and 

^(-^ } )* (E i*) = E :U> Vfc + 1 < j < r - 1, 

and it follows that, 

t: 1 d2c i ) = r(-zwy...r(-z^r (E4q). 

i=l i=l "» 

Using again (20), we have that 

1=1 i=l j 

From this, we deduce that, for 1 < j < k — 1 and a € i£(ci + ... + c&, 1), 

r (- z 0'))*a = r(-z (j ' ) )*a. 

Finally, we conclude that 

i=l i=l U i 

1=1 u 'i 



For the proof of Theorem 3.3, we also need the following result for which a proof is given in 

r r 

Hassairi and Lajmi (2001). For h = hjCj + ^ hij in E, we denote h = hjCj and = 

i=l i<j i=l 

r 

E ^'fc' 1 < i < r - !• 
fc=j+i 

Lemma 4.3. z) For all 1 < j < r — 1, the map from E + into L(E); u i — ► t(u^) is differentiable 
on E + . Its differential in a point u is the map 

h .— ► = 2h {j) U Cj + 4(/i (i) Dc j )(u {i) nc j ) 

The map H : £7 + — > T; it i — > t u is differentiable on E + and 
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H'{u){h) = Y, T(u (0) )T(^ (1) )-^(« (i_1) )^(^)T(u {i+1) )...r(u {r - 1) )P(H) 

3=1 

+ 2T(u W )...T(u^)[2L{u)L(h) - L(uh)} 

where t(u^) = Id. 

Note that for u = e, we get 

r-1 r-1 

K^{e) = 2h^U Cj and H'(e)(h) = £ K}(e) + 2L(fc) = 2[£ h^U Cj + L(h)] 

We are now in position to prove Theorem 3.3. 
Proof of Theorem 3.3. Differentiating (5) with respect to x gives 

a'( x ) = t*g'{tx) + t*g'(t(e - x)). 
Setting x = | and replacing t by 2t in (22) give 



1 /, e 

</ i^t-ej = t w, wnere u = 

Taking t = Id in (22) gives 



g'(te) = t^b, where b = -a'(|) = </(e). 



a '( x ) =g'{ x ) +g'(e-x). 

Inserting this identity back into (22) gives 

i*-y Or) - g\tx) = -[t*"V(e - x) - </(i(e - *))]. 
For any s G (0, 1) and x = te, we have 

g'(sx) = g'(ste) = -t* 1 b = -g'(x). 

Hence for any s € (0, 1), 

t*~ l sf{x) - g'{tx) = -s[t*"V(e - sx) - g'{t{e - ax))]. 
Consequently, on letting s — ► 0, we obtain for any x G fl PI (e — fi) and t € T, that 

«/(ta) = f*~V(*)- 
Now consider u G Then (5) can be rewritten as 

a(x) = g{t u {x)) - g{t u {e - x)). 
12 



To differentiate (25) with respect to u, let us consider the functions 

H : E + — ► L(E); m — > t u and a : L(E) — > E; f — ► /(x), 

so that for u G E + , we can write g(t u (x)) = (g o a o H)(u). 
From Lemma 4.3, we have that 

(goaoH)'{u){h) = g'(t u (x))a'(H(u))H'(u)(h),VheE 

= g'(t u (x))a(H'(u)(h)) 

= < g'(t u (x)),H'(u)(h)(x) > . 

Then, for all he E, 

< g'(t u (x)), H'(u)(h)(x) >=< g'(t u (e - x)), H'(u)(h)(e - x) > . 
Hence by (24) we get 

<tl 1 g'(x),H'{u)(h)(x) >=<t: L g\e-x),H'(u)(h)(e-x) >, 
and for any s G (0, 1), 

<t* u 1 g'(sx),H'(u)(h)(sx)> = <t* u 1 g'(x),H'(u)(h)(x)> 

= < C~V(e - sx),H'(u)(h)(e - sx) > . 

Letting s — ► 0, we get 

< CVC^.^'HW^) >=< CV(e),^'(«)W(e) > • 
Note that if u = e, we have for all x G $7 n (e — 0), 

< <?'(x), H'(e)(h)(x) >=< b, H'(e)(h)(e) >, V/i G £. (26) 

r 

Let 6 = 6jQ + ^ bij be the Peirce decomposition of b defined in (23) with respect to 

i=l i<j 

r 

the Jordan frame (cj)i<j< r . We will show that b = ^6jCj. In order to do so, we consider 

i=i 

x = | + exij with e G]0, ^[ and Xjj G i£y such that ||xy|| = 1. It is easy to see that x = t v (e), where 
r c fc /I 

v = ^2 ~7= + V 2 Cj + ^ £Xi i G ^ + ' ^-^ s hnplies that x G $7 and e — x = | — exij is also 

k=l,k^j * 
in 12. 

Inserting /i = /ijj in (21), we obtain 

H'(e)(h) = 2h ij a Ci . 
13 



With the notations used above and from (1), we have 

H'(e)(h)(e) = hij, 

and 



H'(e)(h)(x) = 1 -h ij+ s<h i3 , Xij >c r 



It follows by (23) that 



< g'(x),H'(e)(h)(x) > = <t* v b, H'{e){h){x) > 

= < b,t- l H\e)(h)(x) > . 



After some standard calculation using (3), we get 



1 



Hence according to (26), we have for all G 

1 

VI -2e 2 

From this we readily deduce that bij = for all 1 < i < j < r. 



< bij, >— < bij, hij > . 



Now let x = t{e), with t E T be in fi. Then from (23), we can write 

g'(x) = f-'b 

i=i 

r 

= t*' 1 (Y,( b k ~ b k+i)(ci + ... + c fc )) 

k=l 

where b r+ \ = 0. 

And using Proposition 4.2, we get 

g\x) = j2(h-b k+1 )(P k (x)y 1 . 
k=i 

By the statement (i) of Proposition 4.1, we have 



di x ) = J2( bk ~ &fe+i)log A fe (x) + c, cGH 
k=i 
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= logA p (x) + c (27) 

where p = b r ). 

Finally, inserting (27) into (5) we get 

a(x) = log A p (x) — log A p (e — x). 

4.2 Proof of Theorem 3.4 

For the resolution of the functional equation (6), we will use second derivatives and the following 
result. 

r 

Proposition 4.4. Let x = ^^XiCi + V^ij be the Peirce decomposition with respect to (cj) of an 

i=l i<j 

element x of E. Then, for q = (qi, q r ) in IR r , 

r r— 1 

J2li X i C i + ^2<lj X ij = [lrP(ci + ••• + C r ) + ^ (q k - q k+1 )P(d + ... + C k )](x). 
i=l i<j fe=l 

Proof. We use induction on the rank r of the algebra. The result is obvious for r = 2, in fact we 
have 

q\X\C\ + q 2 x 2 c 2 + 172^12 = q2{x 1 c 1 + x 2 c 2 + x 12 ) + (qi - q 2 )x 1 c 1 

= [q2P(ci+c 2 ) + (q 1 -q 2 )P(c 1 )](x). 

Suppose the result true for a Jordan algebra with a rank equal to r — 1. We easily verify that 

r r r— 1 

^ftXiq + = <7 r Q^£iCi + ^Zy) + ~ 9r)^iQ + X {qj-q r )xij. 

i=l i<j i=l i<j i=l l<i<j<r—l 

r 

As, for a = ^ ajCj, we have P(a)x = ^ didjXij, then 

i=l i<j 
r 

qiXid + ^ qfj-Xjj = q r P{ci + ... + c r )(x) + 

i=l i<j 



r-2 r-1 



[(q r -i - q r )P(ci + ... + Cr-i) + J2^ k ~ Qk+i)P(ci + ... + c fc )](^ XjCj + ^ 

fe=l i=l l<i<j<f— 1 

This with some standard calculation give the result. 
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x 



Proof of Theorem 3.4. Differentiating (6) with respect to x gives 

a[(x) = t*g'(tx) - t*g'(t(e - x)). (28) 
Differentiating (28) once again with respect to x gives 

a'[( x ) = t*g"(tx)t + t*g"(t(e - x))t. (29) 
Substitute x = \e in the above equation then replace t by 2t to get 

g "(te) =t*~ 1 Br 1 (30) 

where B = |a"(§) = g"{e). 

Observe that taking t = Id in (29) gives 

a'{(x)=g"{x)+g"{e-x). (31) 
Inserting this identity back into (29), we get 

g"(x) - t*g"(tx)t = -[g"(e - x) - t*g"(t(e - x))t]. 

For any s £ (0, 1), change x to sx in the above equation, and use the fact that by (30), we have 

g"(sx) = s- 2 g"(x), 

then on multiplication by s 2 , we deduce that for any s G (0, 1) 

g"(x) - t*g"{tx)t = -s 2 [g"(e - sx) - t*g"{t(e - sx))t}. 

Letting s — ► 0, we obtain, for any x G SI n (e — f2) and t G T 

g"(tx) = f- 1 g"(x)t-\ (32) 
Let us now consider u G E+, then (28) can be rewritten in the form 

ai Or) = t* u g'(t u (x)) - t*J(t u (e - x)). (33) 
In order to differentiate (33) with respect to u, we introduce the functions 

H:E + -^ L(E); u .— t u , vr : L(E) — L(£); / ^ /* and : L(£7) — ► -E 1 ; / 1 — > f(x). 

We have 

(?r o H)'{u){h) = H'(u)(h)*, Vh G 
As, for any u G E 1 .^, we have g'(t u (x)) = (g' o <fio H)(u), then 

fe'^o H)'(u)(h) = g"(t u {x))H\u){h){x). 
16 



We also consider the functions 

* : E + — ► L(E)xE; u i — ► (t* u , g'(t u (x))), where x is fixed, and C : L(E)xE — ► E; (f,z) \ — ► f(z). 
Then one can easily see that 

V'(u)(h) = (H'(u)(hy,g"(t u (x))H>(u)(h)(x)) , 

and it follows that 

(Co *)'(«)(/») = C(*(u))*'(u)(h) 

= C\tl,g\t u {x))){H\u){h)\g"{t u {x))H'{u){h){x)) 

= atlg"(tu(x))H'(u)(h)(x))+C(H\u)(h)\g'(t u (x))) 

= t:g"(t u (x))H'(u)(h)(x)+H'( U )(hng\t u (x))). 

Now differentiating (33) with respect to u gives 

e u g"{t u {x))H'{u){h){x)+H'{u){hy{g'{t u {x))) = tig'' {t u {e-x))H' {u){h){e-x)+H' {u){hf {g' {t u {e-x))). 
If we make u = e, we get 

g"(x)H'(e)(h)(x) + H' (e)(h)* g' (x) = g"(e - x)H'(e)(h)(e - x) + H'(e)(h)*g'(e - x). (34) 
That is 

(g"(x) + g"(e - x))H'(e)(h)(x) + H'(e)(h)*(g'(x) - g'(e - x)) = g"(e - x))H'(e)(h)(e). 
Using (28) and (31), this becomes 

a'{{x)H'{e){h){x) + H' (e){h)* a'^x) = g" (e - x))H'{e)(h)(e). 
On the other hand, we know that 

{a"(e — x) = a"(x) 
a'i{e — x) = —a'i{x). 

Hence 

a'((x)H'(e)(h)(e - x) - H' {e){h)* a'^x) = g" (x)H' \e)(h)(e). 
Given (31), we obtain 

H 1 \e){h)* a'^x) = g"(e - x)H'(e)(h)(e - x) - g" (x)H' (e){h)(x), Vh G E. 

Setting h = e in this equality, we obtain by (21), 
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a' 1 (x)=g"(e-x)(e-x)-g" (x) (x) . (35) 

Therefore, using (35) and the fact that, for all x in n (e — fl), a[(x) = g'(x) — g'(e — x), the 
equality (34) becomes 

/(^'(eX/^-^eX^y^ (36) 
For s 6 (0, 1) change x to sx, to obtain 

g"(x)H'(e)(h)(x)-H'(e)(h)*g"(x)(x) = s[g" (e- sx)H' (e)(h)(e- sx) - H' (e)(h)* g" (e- sx)(e- sx)}. 
Letting s — ► 0, implies that for any x £ £2 l~l (e — f2) and h £ E, 

H'(e)(h)*g"(x)(x) = g"(x)H'(e)(h)(x). (37) 
If x = |, then, using the fact that <?"(§) = 45, we can write 

H'{e){h)*B(e) = BH' (e)(h)(e), V/i G 5. (38) 

Now let 

r 

i=l i<j 

be the Peirce decomposition of B(e) with respect to the Jordan frame (cj)i<j< r . We will show that 

r 

i=i 



In fact, substituting h = /ijQ in (38) gives 

i=i 

r r 

i=i i=i 

In particular, if h = Ci, we get 

B(ci) = QB(e), VI < i < r. 

Consequently 

i-1 

2 

Using (38) again for h = hij, we get 
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B(c l ) = X l c i + ^[J2^+ E A ^ ( 39 ) 

k=l k=i+l 



B{hij) = H'(e)(h)*B(e) 
= 2(hijDci)*B(e) 
= 2(c i ah ij )B{e) 

j-l r 

= < Xij,hij > Ci + Xjhij + 2hij[^2 X kj + X jk ],Vi<j. 

k=l k=j+l 

By symmetry of B, it follows that 

< B(hij),Cj >= - < hij, Xij >= 0, Vi < j. 
This implies that Xij = 0, Vi < j. 
Hence, we have 

{B(ci) = AjCj, VI < i < r. 
B(hij) = Xjhij, VI < i < j < r. 
B(e)=Y j X i c i . 
i=i 

r 

Let x = ^ XiCi + ^2 x ij De the Peirce decomposition of x. By Proposition 4.4, we get 

1=1 i<jr 

r 

i=l i<j 

r-l 

= [A r P( Cl + ... + Cr) + ^(Ajfc - Afc+i)P(Ci + ... + c fc )](x), 

k=l 

that is 

r-l 

5 = A r P(ci + ... + Cr) + ^(A fc - A fc+ i)P(ci + ... + c fc ). 

fc=l 

As for each x in f2, there exists a unique t in the triangular group T such that x = t(e), then 
using (30) and the fact that, for all x in E and for all g € G, P(gx) = gP(x)g*, we can write 

#"(x) = f'Bt" 1 

r-l 

= A r t* _1 P(ci + ... + cv)* -1 + ^(A fe - A fc +i)** _1 P(ci + ... + Cfc )* -1 

k=l 
r-l 

= X r P(t*~\ Cl + ... + cr)) + ^(A fe - A fe+ i)P(t* _1 (ci + ... + Cfc)) 

fc=l 
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Using Proposition 4.2, we can write 

r-l 

g"(x) = X r P(x- 1 ) + £(A* - A fe+1 )P((P fe (x))- 1 ). 
fc=i 

This, invoking Proposition 4.1, implies that 

r-l 

g'(x) = -X r x^ - £(A fc - A jfc+1 )(P fe (x))- 1 + 5 
fc=i 

where 5 £ E. 

And by Proposition 4.1, we have that 

r-l 

g(x) = -A r logdetx - ^(A fc - A fe+ i) log A fc (x)+ < 5,x > +ci 

fc=i 

where ci G IR. 
Hence 

= logA p /(x)+ < (5, x > +ci 

where p' = (p[, ...,p' r ) = (-Ai, -A r ). 

As a^x) = g'{x) — g'(e — x) (See (28)), we obtain 

ai (x) = log A pl (x) + log A pl (e- x) + c 2 

where c 2 G IR. 

Finally, from (6) we get 

a 2 (te) = 2 log A p t(te) + 2ci - c 2 + < 5,te > . 
For any there exists a unique t e T such that y = te. Then 

a 2 (y) = 21ogA p /(y)+ < <5, y > +c 3 , 

where C3 = 2ci — c 2 . 

Acknowledgements We are very grateful to Jacek Wesolowski for discussions on the subject of 
this paper. 
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